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ABSTRACT 



We present results from two observations of the cosmic microwave background 
(CMB) performed from the South Pole during the 1993-1994 austral summer. Each 
observation employed a 3° peak to peak sinusoidal, single difference chop and consisted 
of a 20° X 1° strip on the sky near our SP91 observations. The first observation 
used a receiver which operates in 3 bands between 38 and 45 GHz (Q-band) with a 
FWHM beam which varies from 1° to 1°15. The second observation overlapped the 
first observation and used a receiver which operates in 4 bands between 26 and 36 
GHz (Ka-band) with a FWHM beam which varies from 1°25 to 1°7. The Ka-band 
system has a similar bcamsize and frequency coverage as the system used in our SP91 
results. Significant correlated structure is observed in all bands for each observation. 
The spectrum of the structure is consistent with a CMB spectrum and is formally 
inconsistent with diffuse synchrotron and free-free emission at the 5 a level. The 
amplitude of the structure is inconsistent with 20 K interstellar dust; however, the 
data do not discriminate against flat or inverted spectrum point sources. The root 
mean square amplitude (±1(t) of the combined (Ka+Q) data is AT^ms = 42.0lg^g^ /uK 
for an average window function which has a peak value of 0.97 at £ = 68 and drops to 
g-0.5 peak value at ^ = 36 and i = 106. A band power estimate of the CMB 

power spectrum, Q, gives (^^^^^) = 1.77tJ;|4 ^ 10~^°. The band power estimates 



for the individual Ka and Q-band results are larger than but consistent with the band 
power estimate of the combined Ka-band SP91 results. 




Subject headings: cosmic microwave background — 



cosmology: observations 
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1. Introduction 

Anisotropy measurements of the cosmic microwave background (CMB) are a very effective 
tool for testing and constraining models of cosmic structure formation. With the discovery 
of large angular scale anisotropy by COBE (Smoot et al., 1992), there has been an increased 
interest in characterizing anisotropy on degree angular scales. Although all CMB measurements 
have to overcome a long list of systematic effects and foreground contaminants (Wilkinson, 
1994), they have the potential to constrain many of the global parameters of the Universe and 
thus discriminate among the plethora of cosmic structure formation models. Over the past six 
years, we have travelled to the South Pole three times to perform these degree scale anisotropy 
measurements. The results from our 1988-89 measurements are detailed in Mcinhold and Lubin, 
1991 (SP89a) and Mcinhold et al., 1993 (SP89b) while the results from the 1990-91 measurements 
are detailed in Gaier et al., 1992 (SP91a) and Schuster et al., 1993 (SP91b). The results from 
these measurements and the five balloon-borne Millimeter-wave Anisotropy experiment (MAX) 
are summarized in Lubin, 1994. In order to obtain additional data and frequency coverage, we 
returned to the Amundsen-Scott South Pole Station during the austral summer 1993-94. 



2. Instrumentation 

The SP94 observations used the Advanced Cosmic Microwave Explorer (ACME) as have all 
our previous degree scale anisotropy measurements. ACME is a one meter off axis Gregorian 
telescope which has been described in detail in SP89b. During the observations, the ellipsoidal 
secondary oscillated sinusoidally at 8 Hz with a peak to peak throw of 3° on the sky. The 
receiver signals were phase synchronously demodulated using a "square-wave" lockin amplifier 
and sampled every 0.5 seconds. The beam profile of the telescope can be approximated as a 
Gaussian beam with a 1 cr;, dispersion which varies in frequency as given below. Two different 
total power radiometers were used in these observations. The lower frequency (Ka-band) receiver 
is similar to that described in SP91a and incorporates a very low noise, high electron mobility 
transistor (HEMT) amplifier (Pospieszalski et al., 1990) cooled to 4 K in a "^He dewar. This 
receiver operated at four center frequencies (27.25, 29.75, 32.25, and 34.75 GHz) with 2.5 GHz 
3 dB bandwidths. The band subdivision is used to compensate for gain variations across the 
full band and to obtain spectral information which can be used to discriminate between the 
various astrophysical foregrounds. For the Ka-band system, the beam dispersion is given by 
ah = 0.°70ib0.°04 X ^ JJ^ ^ . The higher frequency receiver (Q-band) is described in Gundersen et 
al., 1994 and also uses a cryogenic HEMT amplifier based on a design developed at the National 
Radio Astronomy Observatory (NRAO). This amplifier was built at UCSB with assistance from 
NRAO and uses an AlInAs/GalnAs/InP HEMT (Pospieszalski et al., 1994) in the first of five 
amplification stages. The Q-band system was multiplexed into 3 equal bands centered at 39.15, 
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41.45, and 43.75 GHz with nominal 3 dB bandwidths of 2.3 GHz, and the beam dispersion is 

given by ai, = 0.°47 ± 0.°04 x (^ ^^'^ ^ • The HEMT amphfiers introduce intrinsic cross correlations 
between the bands which can be characterized by the correlation coefficient between any two 
frequencies. The measured correlations were typically 0.25 and 0.50 for the Ka and Q-band 
systems, respectively, including atmospheric correlations. The radiometers are calibrated to 10% 
absolute accuracy and 3% relative accuracy using a combination of cryogenic cold loads, the sky, 
ambient Eccosorb, and the Moon. The long term stability of the system was checked daily by 
inserting an ambient load "calibrator". These calibrations varied by less than 3% over the time 
scale of an observation and contribute a negligible amount to the final error estimate. 



3. Observations 



Two observations were performed between January 9, 1994 and January 22, 1994 and 
collected 261 hours of data. The first observation used the Q-band receiver and the second 
observation used the Ka-band receiver. These observations consisted of smooth, constant 
declination, constant velocity scans of length 20° on the sky about a center Ucen =45°, Seen =-62°. 
The closest approach to the Sun was 60° on the sky and the closest approach to the plane 
of the galaxy corresponds to b^^ = —40°, = 272°. This is a low foreground emission 
region which overlaps some of the region observed in SP91. Our measurement of the Eta 
Carina region showed that our absolute elevation was one degree lower than we expected. 
The offset has been attributed to sag in the inner frame of the telescope mount and makes 
a direct comparison between these measurements and the SP91 measurements problematic. 
The instantaneous right ascension of the beam for any of the 3 observations can be given by 
a{t) = acen + 20 {S {t) mod 2 — 1/2) / cos Seen + "o sin {2iTVcht) I cos 6cen where S {t) = int {vgct + 1) 
enumerates the scan number, = 1.5° is the sinusoidal chop amplitude, Vch = 8 Hz is the chop 
frequency, and Vgc = 10 mHz is the scan frequency. The instantaneous beam position on the sky 
is then given by if (t) = a{t) cos Seen- Observations of the Moon established the absolute pointing 
at low elevations and this was confirmed with observations of Eta Carina at high elevations. The 
error in absolute pointing is ±0.25° in right ascension and ±0.12° in declination while the error 
in relative pointing is ±0.05° in right ascension and ±0.05° in declination. If a temperature 
at position is compared to a temperature at position n^- at the same declination, then the 
dimensionless window function can be written as 

W,{^^) = B,{at)B,{a^^ Sr^^^^^i^o^K^ " 2r)ao]jo'[(^ " 2r)i^^/2\ cos[(£ - 2r)$^^] (1) 

r=o [2V!(^-r)!] 

where ^>^j = cos~^(n^- n^) is the angular difference (or lag) between the temperature measured 
at bin i with channel k and the temperature measured at bin j with channel I. The beam profile 
function is given by Bi{ab) = exp l)(7^/2] , Hq is the Struve function of index, jo is 
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the 0*'* order spherical Bessel function and = (20/43) (tt/ 180) is the bin size in radians on the 
sky. The indices arc given by i,j = 1 to = 43 bins and k,l = 1 to F = 3 for the Q-band data 
and A;, / = 1 to F = 4 for the Ka-band data. This window function is shown in Figure 2 for the 
different combinations of beamsizes and is a specific example taken from a more general expression 
in White and Srednicki (1994). 



4. Data Reduction and Analysis 

Data were rejected for a number of reasons including poor pointing/chopper performance 
(1.2%), performance of other observations (0.2%), calibration sequences (0.7%), telescope/receiver 
maintenance (9.1%), temperature variations of the cold plate and backend electronics (0.9%), and 
bad weather (12.9%). From a total of 261 hours of data, 196 hours were used in the data analysis. 
The poor weather data were determined in a way similar to SP91b in which the data from a 
single scan were combined into position bins from which an average and one sigma error bar were 
calculated. The of the individual scans was then calculated, and if the probability of exceeding 

for 43 degrees of freedom was less than 0.01 for any channel, then the data from all channels 
for that scan were removed. As a cross check, other weather filters similar to SP91a were also 
implemented with no significant changes in the final data set. As with SP89 and SP91, an offset 
and gradient were removed in time over the time scale of a single scan (100 sees) for each of the 
channels in an observation. The offsets were between 1 and 2 mK depending on the channel. The 
offset and gradient subtraction are taken into account in the analysis by creating a matrix R such 
that = R^J^l where Tj are the F x N temperature means and are the F x {N — 2) projected 
temperatures with a = 1, — 2. The formation of R is discussed in Bunn et al. 1994 and we have 
made R orthogonal such that R^R = I. The coadded means and 1 sigma error bars are shown 
in Figure 1 and show statistically significant, correlated signals for each observation. All quoted 
temperatures have been converted from antenna temperature to thermodynamic temperature and 
have been corrected for atmospheric absorption. In order to determine the origin of the observed 
structure, the data were binned in azimuthal and heliocentric coordinate systems. Various subset 
analyses were performed including dividing the data set into four roughly equal quarters in time 
and dividing the data into four roughly equal quarters depending on the azimuthal position of 
the telescope beam. None of these analyses suggest that the observed structure is anything but 
celestial in origin. 



5. Astrophysical Foregrounds 

There are several astrophysical foregrounds which could contaminate these observations. 
These include diffuse synchrotron and free-free emission from within our own galaxy and 
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extragalactic emission from discrete radio sources. Neither the Sunyaev-Zeldovich effect nor 
diffuse 20 K dust emission is expected to contribute more than a few fiK signal. If we assume 
that the 408 MHz map (Haslam ct al.) is a tracer of diffuse, high galactic latitude synchrotron 
emission, then we calculate the rms differential synchrotron emission to be 1.0 K at 408 MHz for 
these observations. Given the spectral index for diffuse synchrotron is ^ = —2.8 with an antenna 
temperature given by Ta oc i;^, we estimate that the diffuse synchrotron contribution to the 
observed rms to be < 7 /jK in the lowest frequency channel and < 3 fiK in the highest frequency 
channel. The small amount of differential emission that exists in this region of the sky at 408 MHz 
can be correlated with discrete radio sources which have been observed in other source surveys 
at 408 MHz and are identified in the PKSCAT90 database (Wright and Otrupcek, 1990). Unlike 
diffuse synchrotron emission, discrete radio sources and free-free emission cannot be dismissed on 
amplitude arguments alone. For diffuse free-free emission, there are no all-sky surveys which would 
allow a direct estimate of the free-free contribution to the observed structure. Instead, we have 
to rely on estimates based on a 10° x 12° Ha map (Reynolds. 1992) made at a similar angular 
scale to predict an average background free-free brightness temperature that can be expressed as 
T ff = 1 X ^0~^'>Jqh\ CSC b^^ . Reynolds' data suggests that the variations in the Ha intensity may 
be a factor of 2 above the average, such that AT j j = 2T j j . Prom this we calculate the differential 
brightness temperature due to free-free emission (upon the closest approach to the galactic plane) 
to range from ATjj = 30 /iK at the lowest frequency to ATjj = 12 at the highest frequency. 
There have been many discrete source surveys at lower frequencies which arc compiled in the 
PKSCAT90 database. The Parkes-MIT-NRAO (PMN) survey at 4.85 GHz (Wright et al., 1994) 
serves as the most sensitive survey with a flux limit of 30 mJy in our observation region and 
includes all the sources identified in PKSCAT90 for our region. The sensitivity of the telescope 
is 90 /xK/Jy (assuming 100% aperture efficiency), so a 30 mJy flat spectrum source (with a flux 
density S(Jy)(x Tai^'^, P = —2) would produce a 3 /iK signal, which is well below the noise 

of the observations. Since the effective solid angle, Og; of the telescope varies like ilg oc z/^^, a flat 
spectrum point source with a solid angle ilp <C and Tfi^t oc v^"^ would give Ta = f^p-^^^ oc u^, 
while a thermal point source produces an antenna temperature Ta oc i^^. If we make the worst 
case assumption that all the point sources have flat spectra to 45 GHz, then we estimate that 
they would produce a ATj-ms = 30 — 45 /xK. Since the worst case estimates of contamination from 
free-free emission and flat spectrum point sources are comparable to the rms level of the observed 
structure, we cannot dismiss or verify these types of foreground contamination without measuring 
/?. This is addressed in the following likelihood analysis. 



6. Likelihood Analysis 



We use the Bayesian method with a uniform prior (Bond et al. 1991) in the determination 
of the root mean square (rms) amplitude of the data and to make an estimate of the broad-band 



- 7- 



power in the CMB power spectrum (Bond, 1994a, and Steinhardt, 1994). The experimental two 
point correlation function is given by 

, , /Ar(#)Ar(n')\ 1 ^ w 
C{n^ ■ n\) = [—^^if^) = ^EC^e + mWei^^^) (2) 

where (o^mfJ/^') = ^u'^mm'^^ ^ spherical harmonic expansion of the radiation temperature 
given by = Ylitm^imXim{Ti)-, ^£ is the window function (Eq. 1) and Tq = 2.726 ± 0.01 K 

(Mather et al. 1994). The rms amplitude is calculated from 

/ \T\ 2 1 °° 1 

=-E(2^ + l)QT^^ W, = -—Y,W^{^^) (3) 
\ToJrms 47r^ ^pi^ m 

where Npi^ = F"^ x N and Wg is the average window function at zero lag. Following Bond, 1994a, 
the broad band power estimate is given by 

We consider a scale invariant, n = 1, "flat" radiation power spectrum given by 
Ce oc f)^ / (^(^+1)) where the constant of proportionality and fi are determined in the likelihood 
analysis and is the center frequency of channel k normalized to the lowest center frequency of 
the observation(s). The full covariance matrix is given by Mf- = C(n'j^ ■ n^j) + Dfjafaj/TQ where 
D^jafaj is the data covariance matrix which is diagonal in each of the F"^ N x N submatrices. 
This covariance matrix accounts for all spatial and channel-channel correlations as well as the 
error for each bin at each frequency. The offset and gradient subtraction are taken into account 
by creating a projected covariance matrix given by M = RMR'^, where R is defined in Section 

4. The likelihood of the data set is proportional to M exp(— x^/2), where = T'^M'^T 
for F X (N — 2) degrees of freedom. Table 1 lists the resulting rms and band power estimates for 
the individual observations as well as the combined observations. Figure 2 shows the band power 
estimates with /3 = in relation to the COBE band power estimate. 



7. Discussion 



The Q-band broad band power is consistent with standard CDM normalized to COBE and 
has a spectral index, which is consistent with the CMB spectrum. The Ka-band observation is 
also consistent with standard CDM normalized to COBE; however, the spectral index does not 
rule out discrete radio sources and thus does not afford such a straightforward interpretation. 
When the Ka and Q-band data are combined, the broad band power is consistent with the 
standard CDM model normalized to COBE and the spectral index is consistent with the CMB 
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spectrum. For the standard CDM power spectrum, the amphtude of the combined Ka+Q 
observations corresponds to a Qrms-PS = 16.2]'^2 g A*K. The combined Ka+Q data show a marked 
improvement on the 1 a confidence interval for the best fit /3. For the most probable ATrms, 
spectra with j3 < —2 (such as diffuse synchrotron and free-free emission) are formally excluded 
at the 5 a level. The SP94 results are most easily compared to the SP91 results and the results 
of Wollack et al. 1993 (SK93). A reanalysis (Bond, 1994b) of the combined SP91a and SP91b 
observations gives (C(^ ^ = I.OG^q'P x 10^^'^ which compares to the Ka-band SP94 results of 

= 1.17]'^q'42 X lO^^^''. The Bond, SP91 analysis assumes that /? = and that there is no 

channel to channel correlations. The SK93 result gives (c^^)^ = 1.3llo:7 x 10"^°, for /3 = 

The SK93 result is consistent with both the SP94 and the combined SP91 results, although one 

should note the differences between the experiments which are addressed in SK93. 
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The means, covariances and window functions will be made available from rot.ucsb.edu using 
anonymous FTP in the directory [sp94dat] and from the Web site http: / /www. deepspace.ucsb.edu . 
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Tablc 1. Derived Paramoters From Likelihood Analysis (ilo") 



Band 


AT 

^-'^ rms 


/3 


AT 


/ Cil(^+l) \ 
\ 2u Ib 


(1) 


(2) 


(3) 


(//K) 
(4) 


X 10-10 
(5) 


Q-band 
Ka-band 
Ka+Q 


30.8+1^ 
31.9+i\6 


o.2in 


49.3+1^3-^ 
32.011:^/ 
42.0ti:^8« 


2 14+2-00 

1 17+1.33 
^•-•^'-0.42 
■1 77+1.58 
-0.53 



Note. — (cols. [2-3]) The rms amplitude and spectral 
index /3 determined from the Hkelihood analysis using 
a scale invariant, n=l, "flat" radiation power spectrum 
given by Ci oc {f^f'')'^/{i{£ + 1)), where is the center 
frequency of channel k normalized to the lowest frequency 
of the observation; (col. [4]) the rms amplitude for /? = 0; 
(col. [5]) CMB broad-band power spectrum estimate for 
(3 = 0. 
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Fig. 1. — Coadded single difference data in /iK thermodynamic temperature, referenced to the top 
of the atmosphere with ±1 a error bars. The error bars account for the correlated noise between 
the channels. The data in this figure shows that there is significant correlated structure. This data 
set was not explicitly used in the analysis since the spatial correlations between the bins have to 
be taken into account using the full theoretical covariance matrix. 
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Fig. 2. — The top figure shows the band power estimates as given in Table 1 compared to the 
standard CDM radiation power spectrum (Sugiyama, 1994) normahzed to the COBE Q^^l_pg = 
19.9 lb 1.6 fiK (Gorski et al., 1994). The lower figure shows the window function at zero lag for 
the Q and Ka-band receivers. The different traces represent the window functions for the various 
combinations of beamsizes for each receiver. The combined Ka+Q-band window function is not 
shown for clarity sake. 



